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ON CERTAIN TWO-POINT PROPERTIES OF GENERAL FAMILIES 

OF CURVES* 

BY 

JESSE DOUGLAS 

In a paper published in 1903,t R- Straubel proved an important theorem of 
geometrical optics to the following effect. In any isotropic medium, — in 
which the index of refraction is therefore a function of position only, v (x, y , 
z ) , — let a and b be any two points. Imagine an element of area S(Xi to be 
placed at b normal to the (curvilinear) ray ab , and let it subtend at a the solid 
angle 5wi . Reversing the situation, suppose an element dff^ placed likewise 
at a to subtend at b the solid angle 5w2 . Then the relation 

5o-i 60-2 

will subsist. In n dimensions v'^~^ takes the place of v^ . 

Straubel's theorem was made the subject of a paper by T. Levi-Civita,t 
who generalized it to the case of a non-isotropic medium. Levi-Civita's 
fundamental result is as follows. Consider the geodesies of a general (curved) 
n-dimensional variety F„ . Denoting by a and b any two points in F„ , let an 
(n — 1 ) -dimensional element of extent 6<ri placed at b normal to the geodesic 
ab subtend at a the solid angle 6wi . If then an element S(X2 placed likewise 
at a subtends at b the solid angle 6w2 , we shall have 

5wi 6w2 ' 

The object of the present paper is to develop general geometric theorems 
in the above order of ideas for an arbitrary family of curves. In particular, 
the property of Levi-Civita's theorem is shown to hold for a wider class 



* Presented to the Society, April, 1917. See Bulletin of the American 
Mathematical Society, vol. 23 (1916-17), p. 441. 

t Ueber einen allgemeinen Saiz der geometrischen Optik und einige Anwendungen, P h y s i - 
kalische Zeitschrift, vol. 4 (1903), pp. 114-117. 

} Una proprieta de simmetria delle traiettorie dinamiche spiccati da due punti, Atti della 
Reale Accademia dei Lincei, vol. 24, April, 1915. The notation is different 
from that of the present paper. 
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of curve families than the geodesies. Straubel's theorem is discussed as to its 
purely geometrical implications; it is shown, among other things, that the 
Straubel relation is not characteristic of the curve families that can arise as 
the totality of paths of light in an isotropic medium, — that is, of the natural 
families, so called.* 

The method of investigation is characterized by the systematic use of the 
fact that, given a family of curves defined by any system of differential 
equations, the curves infinitesimally near to a given one obey a system of 
linear differential equations — the equations of variation (Poincare) or auxiliary 
equations (Darboux) of those that define the given curve family. 

I. Theory foe two dimensions 

1, Definitions. The functions V {a,b) and W{a,b) . Consider any family 
of CO 2 curves in the plane, as defined by a differential equation of the second 
order, y" = F{x,y ,y'). 

Let a and h be any two points of the plane. Then, under suitable restric- 
tions,! there is a unique curve C of the given family joining a and h . Imagine 
an infinitesimal element hh' of length Scri to be placed at h normal to C, and 
let Ci be the curve of the given family that joins h' to a. Denoting by 6wi 
the infinitesimal angle enclosed at a between Ci and C, the ratio Sai/Scoi is a 
function of the two points a and h ; we call it V {a, h) . 

Similarly, let the element aa' of length 60-2 be normal to C at a , and let C2 
denote the curve that passes through a' and h . Then designating as 6w2 the 
angle at h between C2 and C, the ratio 60-2 /6co2 is appropriately V {h , a) . 
It is to be remarked that in general V {a,h) 9^ V {h , a) . 

Besides the defining equations 

"-(«.'> = £- ^'(»-«'=a^:' 

we introduce also the notation 



W{a,h) = 



V{b,a)- 



Our specific object is the study of the geometric functions V {a,b) ,W{a,b) . 
2. Expression tot V {a,b) . We have under consideration a doubly-infinite 
system of curves, defined by a differential equation of the second order, 

(1) Y" = FiX,Y,Y'). 

To every set of initial values x, y, p for X , Y , Y' there corresponds a 

* After E. Kasner, Differential-geometric aspects of dynamics, Princeton Col- 
loquium Lectures (1912), pp. 34-72, or his papers in these Transactions, vol. 
10 (1909), p. 72, and vol. 11 (1910), p. 121, where several sets of completely characteristic 
properties of natural families are given. 

t See Picard, Traite d' Analyse, vol. 3, pp. 90-99. 
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definite curve; let its equation be 

(2) Y =Y{X;x,y,v). 

Explicitly, we may write 7 as a power series; the first terms are 

(20 Y = y + 'p{X-x)+^F{x,y,v){X-xf+---. 

Let the coordinates of a be {xi,yi) , those oi h , ( X2 , 2/2 ) ; let the slope of 
the curve C be pi at a, and p2 at 6 . Then in (2), C corresponds to the initial 
values xi , 2/1 , pi; also to the initial values a;2 , 2/2 , P2 • Suppose that for these 
values (2) reduces to 

(3) Y = y{X), 

which is then the equation of C . 
The curve Ci will correspond to 

X = xi, y = yi, p = pi + 8pi, 
denoting by the last the slope of Ci at a; the equation of Ci will be therefore 

Y = y{X) + bpi Yp{X; xi, yi,pi) , 

where the powers of dpi beyond the first have been neglected as infinitesimals 
of higher order. The subscript p denotes partial differentiation with respect 
to that variable. 

Since Xi, yi, pi are fixed, Yp is a function of X alone; we introduce for it 
the notation 

rii{X) = Yp{X; Xi,yi, pi) , 

and write the equation of Ci in the form 

(4i) y = 2/(Z) + 6pi'7i(^). 

In the same way, the equation of C2 will be 
(42) Y = y{X) + hp,ri,{X), 

where 

rii{X) = Y^{X; xi, yi, pi) . 

In the figure of the preceding section let hb" denote the element of length 
at b parallel to the y-a,xis and included between the curves C and Ci. Obviously, 

bb" = dpiTjiixi) . 
Hence 

8pi -qiiXi) 



8ai = bb" cos < b' bb" = 



n + piy- 

Also 

Scci = 5 tan~^ pi = :, — , — 2 • 
I + Pi 



292 JESSE DOUGLAS [July 

Therefore 

(5i) ^(^'^)=5^=(l+pi)i/^'?i^^^)- 

To obtain the formula for F ( & , a ) , we have only to interchange system- 
atically the roles of a and b . The result is 

(52) ^^^''')=5^, = -(TT^FF'''^'^^^' 

where the significance of the minus sign is as follows : in forming V (a, b) 
the curve C is viewed in the sense ab ; in forming V (b , a) , in the reverse 
sense ba . 

3. Expression for W{a, b). The important fact concerning vii^) and 
r]iiX) is that they are two independent solutions of a certain linear differ- 
ential equation of the second order.* To see this, return to the equations (1) 
and (2) of the preceding section. (2) represents the general solution of (1); 
this means that we have identically in the four variables X , x, y , p: 

Y"iX; x,y,p) =F[X, Y{X;x,y,v), Y' {X; x,y,p)], 

where the accents denote partial differentiation with respect to X . We may 
differentiate this identity with respect to p , obtaining : 

{Ynv = Fy{Y,) + Fy'{Y'),A 

Under the assumption that all the derivatives involved are continuous, the 
differentiations with respect to X and p are commutative; the above equation 
may therefore be rewritten: 

{Y,)" = Fy{Y,)^^Fy'{Y,y. 

When x,y ,p &ve set equal respectively to xi,yi,pi, Yj, reduces to i?i ( Z ) ; 
when the same letters are equal to x^, yi, pi, Yp reduces to t?2 ( X ) . It 
follows that rii{X) and r]2{X) obey the differential equation 

(6) v"iX) =F,'r^'iX)+F,r,iX), 

where it is understood that in Fy'{X, Y, Y'), Fy{X, Y , 7') we are to 
substitute for Y and Y' from the curve C: 

Y = y{X), Y' = y'{X), 

so that these coefficients are functions of X alone, where X represents the 
abscissa along C . 

The functions t?i, i?2 will be completely defined by giving their values and 
those of their first derivatives for an initial value of X . These are determined 
as follows. Differentiating (2') with respect to p, we have, after writing x, 

*The equation of variation of the differential equation (1). 

t For convenience of notation, small letters are used instead of capitals in denoting the 
partial derivatives of the function F . 
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y, p equal to xi, yi, pi: 

7?i(Z) = (X - xi) + ^ F^x,, y^ , p^) (X - X,)' + ■■■, 



hence 






(7i) 


r?i(a:i) = 0, 


v'li^i) = 1- 


Similarly we find 






(72) 


mi^i) = 0, 


'?2(a;2) = 1. 



Now r]i{X), T?2(X), as two independent solutions of the linear differ- 
ential equation (6), must obey the relation 



T?l(Z) Vi{X) 



Vi(X) V2(X) 



Fy'dX 



where k is an undetermined constant. 
In this write successively 

X = xi , X = X2 , 

and take account of the equations (7i), (72); there results: 

- ■qiixi) = e-^" , 

J Fy'dX 

Dividing the second equation by the first, the unknown k is eliminated: 

??i(a:2) r'Fy'dx 

(8) — r = e-^^i 

172 (a;i) 

The result thus obtained applies directly to the derivation of an expression 
f or ^ ( a , & ) . By definition, 

so that, substituting the expressions (5i), (82) for F ( a , 6 ) , V (b , a) , we have 

wr ,, il+ply "viix2) 

Hence, by equation (8), 

yy{a,o) ^_^p2^3/2e 

This can easily be transformed into 

n(F:'-^MiiL)dx 

(9) W(a,b) = e-^^^y 1+// . 

The function W{a, b) is an intrinsic geometric quantity, being, that is to 
say, independent of the selection of coordinate axes. It is desirable to find 
for it an expression in which only intrinsic elements shall figure. 
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This is attained through the following considerations. We are dealing 
with a family of oo^ curves of which one passes through every point in 
every direction. Let dy/dd denote the rate of change of curvature with 
respect to variation of direction as, remaining at any fixed point, we pass from 
curve to curve of the family. The expression for dy/dd is 

dy _ {l+y'')F,'-dy'F 

36 (l+/)3/2 

and to pass from this to the integrand of (9), we have only to introduce the 
element of arc along C , 

ds = {I + y'^yi''dx. 
Multiplying, we get 

dy 



ids = i^Fy' - -Y-^jnjdx; 



therefore, according to (9), 






(10) W{a,h)=e 

The pdth of integration is the arc ah of the curve C . 

4. The multiplicative property of W(a, b) . From the result (10) certain 
simple corollaries follow directly. 

Consider any three points a, b , c located not in general position but on the 
same curve C of the double infinite family under consideration. We have 

fsy, r^ri rtLri 

W{a,b) = e-^^de""', PF(6, c) = e-^» ae , W{a, c) = e-^- 3e \ 

Multiplying together the first two of these equations and taking account of 
the additive nature of integration, we find 

(11) W{a,b) ■ W{b,c) = W{a,c). 

Along any individual curve C the function W{a,b) is of multiplicative 
character. 

5. Transitivity of the relation V {a, b) = V {b , a) . A second result fol- 
lows as a sub-corollary of (11). We have remarked that in general V{a, b) 
9^ V {b , a) . But suppose, as a special situation, that a,b , c are three points 
on the same curve C , and that 

Via,b) = Vib,a), and V (b, c) = V {c, b) . 

Then it must follow as a consequence that 

Via,c) = Vica); 

for the above relations are equivalent to 

(V{a,b) = l, W{b,c) = l, 
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from which by (11) it follows that 

Wia,c) = 1, 
or that 

V{a,c) = Vica). 

Thus, along any individual curve C , the relation V {a , b) = V (b , a) is a 
transitive one. We remark that this fact is not at all evident from the geo- 
metric definition of V (a, b) . 

6. Curve families with symmetric V {a, b) . The particular doubly in- 
finite systems of curves for which V {a,b) is a symmetric function of the 
two points a and b are found immediately from a consideration of (10). We 
wish to have 

V{a,b)^V{b,a), 
which is equivalent to 

or again, equivalent to 






W(a,b)^l, 



for any two points a and b . For this, the integrand dy/dd must be identically 
zero. Recalling the meaning of dy/dd , we have therefore the 

Theorem. The necessary and sufficient condition that, for a family 
of 00 2 curves in the plane, the corresponding two-point function, V {a,b) , be 
symmetric, is that at any given point the curvature be the same for all the curves 
of the family through that point. 

This means that y may be a function of x,y, but not of y' . Curve families 
with symmetric V {a, b) are therefore characterized by differential equations 
of the form 
(12) y" = il+y''yi'G{x,y). 

7. Straubel's theorem. Geometric interpretation. The theorem of Straubel 
stated at the beginning of this paper follows from applying the formula for 
W {a, b) to the case of a family of light rays. 

Imagine in the plane an isotropic optical medium, in which the index of 
refraction is v{x, y) . Such a medium gives rise to a system of oo ^ curves 
as paths of light, these being determined as the extremals of the calculus of 
variations problem 

v{x, y)ds = lain., 



f' 



or by the equivalent Euler-Lagrange differential equations : 

(13) 2/" = (L, -2/'i.)(l + /) 

where L = log v . Following E. Kasner,* we shall term a family of curves of 
the type (13) a natural family. 

*Princeton Colloquium Lectures, pp. 34-37. 
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Fy' -^-7—;^= - L^-y' Ly= -^Lix,y), 



In applying the formula (9) to natural families, we have 

1 + 2/' 
therefore 

e^^^_v{xi^yi)_ 
^y'^''') - e^^--y^^- v{x.,,y^)' 
that is 

,,,. V{a,h) ^ via) 

^ *^ V{b,a) v{h)' 

which is Straubel's result. 

Straubel's theorem is not purely geometric in so far as it contains the indices 
of refraction v {a) , v {h) . An exclusively geometric form is obtainable by 
considering any three points a, 6 , c of the plane. We shall have 

from which equations the indices of refraction can be eliminated, giving 

W{a,h) ■ W{h,c) = W{a,c). 

Therefore : 

The multiplicative relation (11) which holds in the case of any doubly infinite 
family of curves when a, b , c are on the same curve, holds in the case of a natural 
family for general position of a, b , c . 

This is the purely geometric content of Straubel's theorem. 

8. Converse of Straubel's theorem. Is the Straubel property characteristic 
of natural families? For the answer we return to (14). When can we define 
a function v{x,y) of the points of the plane and have for any two points 

a and b : 

V{a,b) ^ Ha) , 
V(b,a) v{b)- 

By (9), the condition in question is equivalent to 

I \Fv' - 7"T — 72 ]dx = \ogv{xi,yi) -\ogv{xi,y2) -=- L{x, y)\lll\. 
Jxi \ 1 + y / 

For this relation the necessary and sufficient condition is 

„ Sy'F d 

The last is to be regarded as a partial differential equation for F . It can 

be integrated by means of the integrating factor 1/(1 + y' )^ , using which, 

we find 

F L,-y' L, 



(l+/)« (1+/)* 



G{x, y), 
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where G{x, y) is the constant of integration (with respect to y') . Thus we 
obtain for F the expression 

(15) F= iL,-y'L.)(l+y") + G{x,y){l+y'')i, 

and have the theorem : 

The Straubel property is not characteristic of natural families — it holds for the 
more general type defined by (15).* 

9. Theory for an arbitrary surface. In extending the theory of the func- 
tions V and W to a family of oo ^ curves upon a general surface, no new diffi- 
culties present themselves apart from the increased length of calculation. 
The typical feature is that the geodesic curvature yg replaces the y of the 
plane. 

The main result arrived at is the formula for W{a, b): 

W{a,b) = eJ" d8 "^^ • 

We consider particularly the special systems of <» ^ curves for which V (a, b) 
= V {b , a) . Evidently, these are the ones for which dyg/dd = , hence the 

Theorem. The necessary and sufficient condition that a family of <^'^ curves 
upon an arbitrary surface have the property V {a,b) = V (b, a) is that the 
geodesic curvature be the same at any given point for all the curves through that 
point. 

Such families are characterized by equations of the form yg = G{x, y) , 
X and y being any curvilinear coordinates upon the surface. The most inter- 
esting particular case occurs when G{x, y) is identically zero, which condition 
defines the geodesies of the surface. The geodesies of any surface have there- 
fore the property V {a, b) = V {b , a) . This is the content (for n = 2 ) of 
the theorem of Levi-Civita enunciated at the beginning of this paper. Levi- 
Civita's result thus enters as a special case under the theorem stated above. 

II. Theory for ordinary three-dimensional space 

10. Generalities; definitions. In extending the preceding considerations to 
space, we have to consider a family of curves defined by two differential 
equations of the second order: 

Y" = F{X, Y,Z, Y',Z') 

(1) 

Z" = G(Z, Y,Z, Y',Z'). 

Such a family consists of oo* curves, one for each initial point and direction 

* The geometric properties of curve systems of the type (15) have been studied by K. 
Ogura, On the integral curves of ordinary differential equations of the second order of a certain type, 
Tohoku Mathematical Journal, vol. 8 (1915), pp. 93-107. 
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X, y , z, p, q. Let the equations of the curve corresponding to these initial 
values be 

Y = Y{X; X, y, z,p, q) 

(2) 

Z = Z{X; X, y,z,p, q) , 

so that (2) represents the general solution of (1). We can write the equations 

(2) explicitly as power series; the first terms are 

Y = y + p{X-x)+hF{^,y,z,V,(l)iX -xf+ ■■■ 
(2') 

Z= z+ q{X -x)+\G{x,y,z,p,q) {X-xf+ ■■■ . 

If we select any two points of space a and h , there will be, under the proper 
restrictions,* a unique curve C of the family (1) passing through these points. 
Let an element of area 5cri be situated at h normal to C; the curves that join 
the boundary of 5cri to a form at a a curvilinear cone of infinitesimal solid 
angle 5wi . Then we define : 

A similar definition determines V {h , a) , and as in the theory for the plane 
we introduce also the function 

In order to derive formulas for V and W we refer the neighborhood of a to a 
rectangular axis system ^ , ij , f of which the f-axis is the tangent to C at a . 
Similarly, the neighborhood of b is referred to an axis system u, v, w where w 
is the tangent to C at 6 .f We designate the plane ^ , 77 as tti , the plane u , 
« as 7r2 . 

The curves of the given family as they stream away from a and impinge 
upon the plane ttz set up a correspondence between the directions through a 
and the points of X2, — a correspondence which for directions infinitely near 
to f and therefore for points infinitely near to b is linear (affine). Now V {a,b) 
is equal to the determinant of this linear correspondence. For V (a, b) may 
be interpreted as the ratio of two areas correlated with one another affinely, 
and under an affine transformation all areas are multiplied by the determinant. 

IL The functions 77 (Z), r(Z) and their linear differential equations. 
Let the Cartesian coordinates of a and b be respectively (xi, 2/1 , Zi ) , 
{x-i, 2/2 , Zi) . Suppose the equations of C to be 

(3) Y = y{X), Z = z{X); 



* See Picard, Traite d' Analyse, vol. 3, pp. 90-99. 

t In the choice of each axis system there is, of course, one degree of freedom, but this cir- 
cumstance is indifferent for what follows. 
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then the direction of C at a is defined by the parameters 

Vi = y' {xi), qi = 2'(xi), 
and at h by the parameters 

Vi = y' i^i), qi = z' {xi). 

The equations (3) of C will result upon the substitution of X\, y\, 21 , pi , gi 
iov x,y , z, -p , qm the equations (2) ; also upon the substitution of X2, 1/2, Z2, 
fi , qi for the same letters. 

Let Ci be a curve, of the family under consideration, passing through a in a 
direction jpi + 5pi, gi + bqi infinitely near to that of C The equations of Ci 
result by writing in (2) : 

a; = xi , y = yi, z = zi, p = Pi + Spi, g = gi + 5gi . 

When we then expand according to powers of dpi , dqi we may omit the powers 
beyond the first as negligible in the limit, and find for the equations of Ci : 

Y = y(X) + 8pir,i{X) + dqir,2{X) 

(3i) 

Z= 2(Z) + 5piri(Z) + 5gir2(Z), 

where we have written 

?7i(Z) = Yp{X; xi, 2/1, zi, pi, gi) 

(4i) 

ti{X) = Zj,{X; xi, yi,zi, pi, gi), 

172 (X) = Yg{X; xi, yi,zi, pi, gi) 

(42) 

tiiX) = Zg{X; xi, yi, Zi, pi, gi). 

In the same way the equations of the curve C2 passing through 6 in a direc- 
tion Pi + Spi , qi + Sg2 will be 

F = 2/(Z) + 5p2'73(Z) + Sg2i74(Z) 

Z= z{X) + bpiU{X) + bqiU{X), 

where 

'nzi.X) = Yp{X; xi, yi, zi,pi, qi) 

(43) 

h{X) = Zp{X; xi, yi, Zi, pi, qi) , 

i74(Z) = Y^{X; xi, yi, zi,pi, qi) 

(44) 

U{X) = Z^{X; xi, yi, zi, pi, qi) . 

The four systems 77, f just defined are independent solutions of a certain 
system of two linear differential equations of the second order.* For since (2) 
represents the solution of the equations (1), there will result upon substitution 

* The equations of variation of the differential equations (1). 

Trans. Am. Math. Soc. 30 
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of (2) in (1) two identities in X ,x, y , z,p , q. Differentiating these identities 
with respect to ^5 , we obtain : 

( Y")p = FyY, + F. Z^ + Fy'(Y')^ + F.'iZ')^ 

{Z")p = Gy Yp -\- GzZp + Gy'{Z')p + Gz'{Z')p. 

Assuming the derivatives involved to be continuous, the differentiations with 
respect to X and p are commutative, and the above equations may be re- 
written : 

{Yp)" = Fy{Y,) + F,{Zp) + Fy'{Yp)' + F^'iZp)' 

{Zp)" = Gy{Yp) + GAZp) -{■ Gy-iYp)' + G^Zp)' . 

For x,y,z,p, g equal to a;i, yi,zi, pi,qi, Yp, Zp reduce to 771 (Z), ^i{X). 
It follows that rji , f 1 obey the linear differential equations 

V'(Z) = Fy' v' (X) + F.' nX) + FyviX) + F. r (X) 

^^^ nX) = Gy'v'{X) + G,'t'{X) + Gyr,{X) + G.nX), 

where it is understood that, in the coefficients Fy' {X , Y , Z , Y' , Z') , etc., 
the substitutions 

Y = y{X), Z = z(X), Y' = y'{X), Z' = z' {X) 

are to be made, so that these coefficients are functions only of X , the abscissa 
along C . 

In similar manner the other systems 7? , f can be shown to obey the equations 
(5). Each ?7, r is therefore determined if initial values are given for 17, f, 
7]' , f ' . For 171 , f 1 such values are found by differentiating (2') with respect 
to -p , and then replacing x , y , z , p , qhj xi , yx , zi , pi , qi \ this gives : 

7ji(Z) = {X - xi) + \ Fp{xi , yi , Z). , fi , qi ,) {X - Xif + ■■■ 

ri(Z) = I Gp{xi, yi,zi, pi, gi,) (Z - .ti)^ + • • ■; 

hence 

(61) 77i(a;i)=0, ri(.'J;i)=0, r,[(xi)=l, ri(^i)=0. 

Similarly we find 

(62) ViiXl)=Q, f2(.Tl)=0, v'2iXl)=0, r2(.T;i) = l> 

(63) Vzix2)==0, hix2)=0, rj',{x2) = l, tUx2)=0, 

(64) 5?4(a;2) = 0, f4(a;2)=0, ri'i(x2)=0, fUa;2) = l. 

The fact important for the present developments is that the four systems 
?j, f , as a fundamental system of independent solutions of the equations (5), 

* For convenience of notation, small letters are used instead of capitals in denoting the 
partial derivatives of the functions F and G . 
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must obey the relation 
(7) 



Vx{X) v^{X) vziX) r),{X) 

ri(z) u{X) u{x) u{x) 

r,[(X) V2{X) vUX) v'UX) 

r;(Z) t2{X) nix) mx) 



r{Fy' + G.')dX 



where h is an undetermined constant. In (7), put X = xi; then by reference 
to (6i), (62), 



(81) 



iizixi) ■rnixi) 



= e-^* 



r\Fy' + G,')dx 

^ k 



Similarly, by putting X = Xi and referring to (63), (64), we have 

f\Fy' + G,')dx 



2) 



ri(a;2) r2(a;2) 



= (jJ k 



On dividing (81) into (82), A; is eliminated; we find 



(9) 







fsCari) 


Vi{xi) 

UiXl) 



= e-^a^l 



This result will be directly useful in deriving an expression for W{a, b) . 

12. Expression for V {a, b) . Recurring to § 10, we next deduce an expres- 
sion for the linear correspondence there described. 

The curve Ci which passes through a in the direction pi + 8pi, gi + dqi 
will pierce the plane X2 in a point 61 of coordinates x^ + 8x2 , 2/2 + Sy^ , 22 + 622 • 
These coordinates must satisfy the equations (3i) of Ci; therefore 

2/2 + 52/2 = y (xi + 5a;2) + Spi rn (x2 + 8x2) + 5f/i 772 (x^ + 8x2), 

22 + 522 = 2 (a;2 + 5a;2) + 5pi tii^i + 8x2) + 8qi r2 (a;2 + 5x2) , 

whence, omitting infinitesimals of higher order than the first, 

8y2 = P2 8x2 + 5pi 7?i (a;2) + 5gi 1)2 {X2) 

522 = g2 8x2 + 5pi fi (3:2) + 5gi r2 (3:2) . 

Let the direction cosines of the u , v axes at b be 

u: a, b, c 

v: I, m, n; 



(10) 



those of the w-axis are 



iv: 



P2 



?2 



^l+pl + qV ^1 + pi + qV ^l+vl + ql 
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Since u, v, ic are a trirectangular system, we have the relations 

(11) 1 + p2m+ Qin = 

bn — cm = 



^i + vl + qf 

At the same time we may turn our attention to the ? , r? , T system at a , 
where, writing for the direction cosines 

of ?: a, iS, 7, 

of 7? : \, IX, V, 

Pi qi 



on-. 



^l + pl + qV Vr+pTrg!' yll + pl + qV 
we note the relations 

(12) 'K + piix + qiv = 0, 

pv - yix = -1======= . 

The coordinates u , v oi hi in the plane ir^ are connected with the space 
coordinates bx2 , hy^ , hzi as follows : 

u = ahx2 + hhyi + cbz^ 

V = I8xi + m8y2 + ndz^ . 
Hence by (10), 

u= {a + p2b + q2c)8x2+ [brjiiXi) + cti{x2)]dpi + [bri2{x2) + c^2(,X2)]Sqi 

V = (I + P2 m + q2 n) 5x2 + [ mrji ( X2 ) + Jifi ( a;2 ) ] 8pi 

+ [mr]2{x2) + nt2ix2)] 5gi; 

but by (11) the coefficients of 5.^2 are zero, so that the equations connecting 
u , V with 5pi , Sqi are 

u = [brii(x2) + cti{x2)]Spi + [bri2(x2) + c^2(a:2)]5qi 
(13) 

V = [mviixi) + n^i{x2)]8pi + [m7]2(x2) + nf 2 ( a;2 ) ] 5gi . 

For temporary expediency we call the coefficients in (13) A, B, C, D, 
writing these equations as 

u = Adpi + B8qi 
(13') 

1) = C8pi + D8qi . 
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The next step is to pass from 8pi , 8qi as coordinates for directions through 
a to ^, ri . ^ is the cosine of the angle between 1 : pi + Spi : qi + 3gi and 
a, j3 , y , hence 



^ = 



a + Pipi + dpi) + jiqi+ 8qi) ^8pi + y8qi 



Vl + p? + ql 



vr 



pi + q'l -^i + pi + qi 

on account of the first equation (12). In the same way, 

\ + fji(pi + 8pi) + v(qi + 8qi) fxdpi + vdqi 



Vr+ pi + q! 

by the second of the equations (12). 
Solving for 8pi , 8qi , we have 



Vl + pf + q! 



8pt = il+pl + ql)(vi-yv), 
8qi = (l+p? + g?)(-Mf+/3^), 

wherein we have used the third equation (12). 

Carrying these values for 8pi, 8qi into (13'), there is obtained for the linear 
correspondence between ^, -q and u, v: 



(14) 



u 



(1 + p? + q\)[{vA -fiB)^- {yA - pB)r,}, 



v={l+pl + qi){{pC - tJ.D) ^ - {yC - fiD)rj]. 



V /x 

y /3 



A B 
C D 



We have said that V{a, h) is the determinant of this correspondence. So 
V{a,h) = {l+p\ + q\y 
Introducing the expansions of ^, 5, C , D as in (13), we obtain 
V{a,h) = {l+p\ + q\f 



V n 

y ^ 



b c 
m n 






and then, by reference to (11) and (12), 

n + pi + qiy" 



(15i) 



Via,b) = 



(l+pi + ql) 



2^1/2 






From (15i) F( 6, a) is at once deducible by the systematic interchange of 
a and b: 



(152) 



Vib,a) = 



(l+pl + qi) 



3/2 



(l+P? + ?f) 



2\l/2 



^sixi) f4(a;i) 



13. Expression for W(a, b) . The preceding formulas show that the deter- 
mination of V {a, b) requires in general the solution of a system of two linear 
differential equations of the second order. But, as in the two-dimensional 
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theory, it appears that W(a, b) , on the contrary, can be expressed by a 
Quadrature. 



quadrature. 
For, by division of (15i) by (152), 



Wia,b) = 



V{a, b) 
V{b,a) 



(l+pl + qlf 


1?l(.T2) 

ri(a;2) 


172 (3:2) 


{1+pl + qiy- 


rsCari) 


riiixi) 



and it has been shown at the end of § 11 that 



Vl ( X2 ) 1?2 ( a.*2 ) i 

tiixi) t2{x2)' r\Fy' + G/)dx 

'73(3:1) Vi(xi) 



Therefore 



(l+v f + q!)' lyPy' + G^') 



+Pl + ql) 



(16) 



Wia,b) = 
An equivalent form is : 



dx 



W{a,b) 



y'F + zG' 



>) 



dx 



14. The invariance of ( 57i/5^i ) + ( dy^/dd^ ) . It is desirable to obtain an 
expression of W(a,b) in geometrically intrinsic terms. For this, we must 
find a quantity to play the role of the dy/dd of the two-dimensional theory. 
Such a quantity is furnished by the following considerations. 

Given any family of =0^ curves in space, 

y" = Fix,y,z,y',z') 

z" = G{x, y, z, y' ,z'), 

each point and direction (lineal element) define a curve. Let 0, Ot he any 
point and direction of space, determining the curve C . 

Construct a rectangular axis system Ot, Ou, Ov of which Ot is one of the 
axes. In this, it is to be emphasized, there is one degree of arbitrariness, for 
the axis system may be rotated freely about Ot . 

The u and t axes define a pencil of directions, and accordingly a pencil of 
curves. These form a surface Si tangent to the plane ut . In the same way, 
the V axis defines with the t axis a pencil of curves, which form a surface S2 
tangent to the plane vt . 

The surfaces Si and S2 intersect in the curve C . Introducing here the idea 
of geodesic curvature, we direct attention to the fact that the geodesic curva- 
ture of a curve upon a surface is relative to both surface and curve. Let the 
geodesic curvature at of C be 71 upon Si , and 72 upon S2 • 
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Regard Si . It is covered with a pencil of curves issuing from . Let 
dyi/ddi denote the rate at which the geodesic curvature changes in this pencil 
as we pass from C to the infinitesimally adjacent curves. And let dyt/ddi 
denote the corresponding quantity for the surface S^ . 

Then we state the 

Theorem. The sum {dyi/ddi) + {dyi/ddi) is independent of the one 
degree of freedom involved in the choice of the triredangular system t , u, v. 

To prove this, suppose the curve family referred to t, u, v as x, y , z axes; 
then, taking advantage of the fact that geodesic curvature is equal to tangential 
ciu-vature, we have 



371 _, dj: 



2 



- = F,' -^, = G. 



ddi " ' de 



2 



If the t, u, V system is rotated about t through an angle (p, we find that, 
relative to the new trirectangular system, 

-^^ = cos^ <pFy' + sin <p cos <p(F/ + Gy') + sin^ (l>G/ , 
dui 

dy2 

-rr- = sin^ <pFy' — sin <p cos <p{F/ + Gy') + cos- <pGz' ■ 

Adding, the angle (p disappears: 

whence the invariantive property of ( dji/ddi ) + ( 572/^^2 ) as stated. 
With this result established, the formula (16) can be given the form: 

„;, ^,„,,).^ilJ^ = /(g+t)^ 

V(o,a) 

15. The case of symmetric V (a, b) . The families of 00* curves in space 
for which F(a, 6) is symmetric in its two arguments are found immediately 
from considering the formulas (16), (17). 

Theorem. The curve families for which 

V{a,b) = V{b,a) 
are those for which 

^1^^1 = 0* 
ddi "^ 3^2 ~ 



* A more graphical form of this condition is as follows: We can construct an indicatrix to 
show the variation of curvature; namely, 0, Ot being any point and direction, by laying off 
on every ray Ou perpendicular to Ot a length proportional to {dy/d6 )~i/^ relative to the pencil 
defined by Ou, Ot. The indicatrix is in every case a conic with as center. When V {a, b) 
is symmetric every indicatrix is an equilateral hyperbola, and conversely. 



F,' + G.'-i^, X^;, -0. 
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How is the property V (a,b) = V (b , a) to be recognized in a given 
family of oo^ curves from the analytic form of the differential equations? 
According to (16) the necessary and sufficient condition for the property in 
question is 

y'F + z'G 

1 + 2/'' + z''- 
This equation is transformable into 

±\ I ] + ±\ ^ 1 = 

Hence the result : 

The families of <x^ curves that have the property V (a, b) ^ V (b , a) are 
characterized by differential equations of the form 

y" = (l+/ + 2'^)20,, 

z" = {l + y'' + z''yQ2, 
where Qi , ^2 are any two functions of x, y , z, y' , z' that obey the relation 

^^^) dy'+dz'-^- 

The general solution of (18) is 

fii = d^''^(^' V' ^> y' ' 2')' 

^2 = - q:j4>{x, y,z,y', z'), 

where is any function of its five arguments. The type of curve system here 
under consideration depends therefore on one arbitrary function of five argu- 
ments, whereas the general system of oo < curves in space depends on two such 
functions. 

16. Straubel's theorem; its converse. Straubel's theorem for three-dimen- 
sional space results from applying the formula (16) for W(a, b) to the case 
of a natural family of curves. Such a family, consisting of the paths of light 
in an isotropic optical medium, is determined as the system of extremals of 

I V {x , y , z) ds , V being the index of refraction, — or by the equivalent differ- 
ential equations 

y" = {Ly-y'L.){l+y'' + z''), 

(19) 

2" = {U-z'L,){l + y'' + z''), 

where L = log v . 



* 



' Cf. E. Kasner, loc. cit. 
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Applying (16) to the equations (19), we find 

^ ' V{b,a)~ p-'ib)- 

Consider the converse question. In connection with what systems of oo* 
curves can a function of position v be defined so that the relation (20) shall 
subsist? By use of (16), it turns out that this is possible whenever the differ- 
ential equations are of the form 

y" = {Ly - y'i,) (1 + / + z'') + (1 + / + /)2 Qi, 
(21) 

Z" = (i, - 2' i,) (1 + / + 3'') + (!+/+ Z'')2 %, 

where fli , Q,^ are any two functions oi x, y , z, y' , z' subject to the condition 

dQ.1 dQ,2 _ 

The Straubel property is therefore not characteristic of natural families. 

17. Flat (Euclidean) n-dimensional space. The preceding theorems admit 
of easy generalization to a flat n-space. Using as notation for the coordinates 
^y yi> 1/2) • • ■ ) Vn-i , a system of n — 1 differential equations of the second 
order, 

(22) y'; = Fi{x,y,-, y'j) {i,j = l,2,---,n-l) 

will define a system of oo^"-^ curves. Such a system will have the property 
V {a, h) = V (b , a) when the differential equations are of the form 

(23) 2/'/ = (l+2/;'+---+2/«-i)<"+""fii (i = l,2,.-.,n-l), 
where the Q's are any functions of the x , y, , y'j that obey the relation 

(24) ^^ + |"-^+...+f- = 0. 
02/1 02/2 oy„-i 

The Straubel property, which for ra-space is 

V{a,b) ;^"-i(a) 
V{h,a) ~ v^-^{b)' 

is found to be characteristic of the families defined by differential equations 
of the form : 

2/7 = ( iyj - 2/i ix) ( 1 + 2/1 + • • • + 2/1-1 ) 

(25) 2 2 

+ (l+2/i + ••• +2/n-i)^"+"''fii 

a = 1,2, ■■■,n-l), 
where the fl's are subject to the condition (24). 
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We note that the general solution of (24) is obtained by taking any ( n — 2 ) 
functions 

(t>l, <^2 , • • • , 071-2 

of the arguments x , y, , y'j and setting the O's equal to the jacobians:* 

d {4)1, (j)2, ■ ■ • , (t>n-2 ) ^ 5 ( <^1 , 02 , ■ • • , 4>n-2 ) 

^^1 — a / ' ' ' — x ) ih = — ^ , ' > } — ;- , etc. 

0(2/2, 2/3, ••• , 2/,.-i) d{yi, 7/3, ••• , 2/,.-i)' 

The families with symmetric V (a, b) depend therefore on one less arbitrary 
function of the 2n — 1 arguments x, y,, y'j than the general family of 



III. Theoby for a general (curved, non-euclidean) variety 

OF n DIMENSIONS 

18. The invariance of '2{dyi/d9i). Formula for Wia,b). We extend 
our developments, finally, to a general variety F„ , one in which the squared 
element of length is defined by the general quadratic form, 



ds^ = X dadxidx 



k- 

The first step is to generalize the theorem of part II concerning the invariance 
of (971/5^1) + {6^2/ ddi). We are considering a system of oo^™-^ curves 
in Vn, one for each point and direction. Let the point and direction Ot 
define the curve C . Then let a system of n rectangular directions be con- 
structed of which Ot is one: Ot, Oh, ■ ■ ■ , Ot„-i. Each Oti determines with 
Ot a pencil of directions, therefore a pencil of curves, which form a two-spread 
V2''. There arise in this way n — 1 Vi's, all intersecting in the curve C . 
On each V'i^ , C has at a definite geodesic curvature yi , which changes at a 
definite rate dji/ddi as we pass from C to the infinitesimally adjacent curves 
upon V'^\ 

Now, in the construction of the rectangular system Oti it is seen that there 
exist |(n — l)(n — 2) degrees of freedom. Hence the significance of the 

Theorem. The sum 

is independent of the ^ (n — 1 ) ( n — 2 ) degrees of freedom involved in the con- 
struction of the rectangular system Ot, Oti, • • • , Otn-i . 
Upon this result we can then establish the theorem: 

(1) W(a,b) = e'^'^-A'^O' . 

From the form of the last equation there follow for a general space the same 
* See De La Valine Poussin, Cours d' Analyse, vol. 2, p. 302. 
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corollaries as in §§ 4, 5 of the theory for the plane. It follows also that the 
necessary and sufScient condition for 

V{a,b) = V{b,a) 
is 

19. Geodesies. Levi-Civita's theorem. Consider the most important ex- 
ample of co2ra-2 curves in F„ — the geodesies. A geodesic in F„ is a fortiori a 
geodesic upon any sub-variety Vm upon which it may lie. Recalling the 
meaning of yt as defined in the preceding section, we see that for the geodesies 
we have, separately and identically, 

7i = , 72 = , • • • , yn-i = , 

and hence, separately and identically, 

^ = ^ = ••■ ^^^ = 

From these equations there follows, according to (2), 

Levi-Civita's Theorem.* The geodesies of any n-dimensional variety have 
the property V {a, b) = V {b , a) . 

Levi-Civita's result thus falls as a special case under the theorem stated 
above in connection with the equation (2). 

20. Conformal representation. Straubel's theorem and its generalization. 
In conclusion, we consider Straubel's theorem for a general variety F„. 
Taking the Straubel relation 

V{a,b) v^-^a) 



(3) 



V{b,a) v'^-'ib)' 



we inquire as to the most general family of w^" ^ curves in F„ in connection 
with which a point-function v can be defined so that the stated relation shall 
hold. 

To this end, conformal representation furnishes the easiest method. Let 
the manifold F„ , of which the fundamental form is 

ds^ = X) ^»fc dxi dxk , 
be represented upon a second manifold F„' of which the fundamental form is 

dsi = j'^ X) «»■* (J^i dxk 

so that points with the same coordinates correspond. Since dsi = vds , the 
mapping so defined is conformal with modulus of conformality v . 
* See the introduction. 
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The behavior of V (a, b) under the conformal transformation is seen to be 
as follows: 

(4) V{auh) =v-'{b)V(a,b), 

(5) Vib^,a^) =y--^a)V(b,a). 
By division, 

(a^ Vjauh) _ v^^Hb) V (a , b) 

^^ V{b^,a^) p--^{a) V (b , a)' 

Hence if in F„ we have the condition (3), then in Vl we shall have 

(7) F(ai,6i) = F(6i,ai), 

and conversely. 

The relation (7) holds in particular for the geodesies of V'„ (Levi-Civita's 
theorem). These are the extremals of 



f' 



dsi = min. ; 

the image curves in F„ are therefore the extremals of 

vds = min.; 



S' 



they are, in other words, a natural family. Thus all natural families have the 
Straubel property. 

But the property (7) belongs not only to the geodesies of V^, — it holds for 
every system of oo^"-^ curves that satisfy the condition I!, (dji/ddi) = 0. 
Hence the Straubel property is not restricted to natural families, that is to the 
conformal images of geodesies, but is true of the more inclusive type that 
results by conformal mapping from any family with the symmetric property 
V (a, b) = F(&,a),orits equivalent S ( dji/ddi ) = . 
Columbia University 



